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Predictions for total inelastic cross-sections for photon induced processes are dis- 
cussed in the context of the QCD-inspired minijet model. Large theoretical uncer- 
tainties exist, some of them related to the parton distributions of hadrons in impact 
parameter space. A model for such distribution is presented, based on soft gluon 
summation. This model incorporates (the salient features of distributions obtained 
from) the intrinsic transverse momentum behaviour of hadrons. Under the assump- 
tion that the intrinsic behaviour is dominated by soft gluon emission stimulated by 
the scattering process, the b-spectrum becomes softer and softer as the scattering 
energy increases. In minijet models for the inclusive cross-sections, this will counter 
the increase from Oj e t . 

The impact of parton scattering on the rise of inclusive cross-sections with energy was 
suggested by Cline and Halzen Q, after such rise was first observed in proton-proton 
collisions at ISR . Minijet models were put forward to describe quantitatively the further 
rise at higher energy 0, [|, |j] and eikonal minijet models ||, || |7|, || were subsequently 
developed to include an increasing number of partonic collisions in QCD resulting from 
the rapid rise in gluon densities. Recent measurements of photo- and hadro-production 
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total cross-sections P, [L0| in energy regions where QCD processes dominate, confirmed 



the rising trend and have been confronted with theoretical predictions obtaining varying 



degrees of success 11 



The simplest mini-jet model || was written as 

&inel{s) = Vsoft + Vjet{Ptmin, s) (1) 

with the rise controlled only by an energy dependent parameter, namely ptmin, which 
regulated the Rutherford scattering divergence in the QCD jet-cross-section. The lack of 
unitarity of this model was amended in the eikonalized mini-jet model, where inelastic 
hadron-hadron cross-sections are written as 

o md = J d 2 b[l-e- n ^}. (2) 
Here the average number of collisions at impact parameter b is given by 

n{b,s) = A{b)a{s) (3) 

with fixed ptmm- In this version, the excessive QCD rise, controlled by ptmin(s) in the non- 
unitarized models, was softened through eikonalization and introduction of the overlap 
function A(b), which thus became a key ingredient of all models with a QCD component. 
This function, which describes matter distribution in impact parameter space, in most 
applications [H] has been assumed to be the Fourier transform of the product of hadronic 



form factors of the colliding particles. In other models |T2|, a gaussian shape has been 
preferred, thus relating A (b) to the intrinsic transverse momentum distribution of partons 
in the colliding hadrons. In either case, detailed information on A(b) relies on parameters to 
be determined case by case. However, while direct measurements of the EM form factors 
are available for nucleons and pseudoscalar mesons, experimental information regarding 
photons or other hadrons such as vector mesons is lacking. The same observation applies 
to the intrinsic-p 4 interpretation for the spatial distribution of partons in vector mesons 
or photons. Thus the extension of this model to photonic cross-sections unveils one of 
its main drawbacks, viz. lack of a fundamental description of parton b-distributions. To 
reduce the uncertainties in the QCD description of the rise of the inelastic cross-section, 
and allow this model to graduate to QCD respectability, it is mandatory to arrive at a 
QCD description of the overlap function A(b). 

For the case of photonic processes, there are further uncertainties related to the hadronic 
behaviour of photons. The model has been adapted to photonic processes by writing the 
inelastic cross-section as 



,nel = phad f ^ _ g-n(M)] ( 4 ) 



where P a 1 b ad gives the probability that both colliding particles a, b be in a hadronic state ]T3"| 
and n(b, s) = n so f t (b, s) +nh ar d(b, s). Here n so f t {b, s) contains the non-perturbative part of 
the cross-section from which the factor P^ has already been factored out and the hard, 
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QCD contribution to the average number of collisions at a given impact parameter b is 
given by 

n har d(b, s) = Aabitf—^aif (5) 
a 3 a l is the hard part of the cross-section. We then have 

rjhad Tjhad. phad / phad\2. rjhad -r 

IP 7 ' 77 ^ 7 ' ' PP 



The predictions of the eikonalised mini-jet model |L4[] for photoproduction processes 



therefore depend on 1) the assumption of one or more eikonals 2) the hard jet cross-section 
a at = Iptmin Iffidpt "which in turn depends on the minimum p t above which one can expect 
perturbative QCD to hold, viz. ptmin-, and the parton densities in the colliding particles a 
and b, 3) the soft cross-section a s °J l 4) the overlap function A ab (b), usually written as 

A ab (b) = J^y 2 J d 2 qF a (q)F b (q)e^ b (6) 

where F is the Fourier transform of the b-distribution of partons in the colliding particles 
and 5) last, but not the least, P^. 

To study the parameter dependence of this model, one can restrict attention to a single 
eikonal : more eikonal terms although improving the fits, de facto introduce new sets of 
parameters and very much reduce the predictivity of the model. The hard jet cross-sections 
have been evaluated in LO perturbative QCD using two different photonic parton densities 
DG ][H| and GRV (TBI. The dependence of a J a l on ptmin for DG densities is given in Ref. 



11]. Clearly this dependence is strongly correlated with the parton densities used. Here 



we shall only show the results for eikonalised mini-jet cross-sections using GRV densities. 
Further, for the purposes of this note, we try to estimate cr]oft fr° m a 1o}t which in turn 
is determined by a fit to the photoproduction data. For 77 collisions, we use the Quark 
Parton Model suggestion cr]^ t = §0"Jo/ t - 

In the original use of the eikonal model, the overlap function A ab (b) of eq. ^| was 
obtained using for F the electromagnetic form factors and thus, for photons, a number of 
authors ]17| [l!| have assumed for F the pion pole expression, on the basis of Vector Meson 
Dominance (VMD). As mentioned, another possibility is that the b-space distribution of 
partons in the photon is the Fourier transform of their intrinsic transverse momentum 
distributions. For protons this has been assumed to correspond to a gaussian shape. 



For photons, the perturbative part [TJ| of the intrinsic transverse momentum has been 



suggested to correspond to the functional expression 

(7) 



dN-y 



dk t k t + k% 

Recently this expression was verified by the ZEUS EO] Collaboration, with k Q = 0.66 ± 
0.22 GeV. It is interesting to notice that for photonic collisions the overlap function will 
have the same analytic expression for both ansatze for F: the VMD inspired pion form 



factor or the intrinsic transverse momentum; the only difference being that the former 
corresponds to a fixed value of ko = 0.735 GeV whereas the latter allows to vary the value 
of the parameter ko- Thus both possibilities can be easily studied by simply changing ko 
appropriately. The overlap function, which for proton-proton collisions would be given by 

as proposed by L.Durand et al. || in the first eikonal mini-jet model for proton-proton 
collisions, for 7J9 collision would become 



.4 



FF 



> 2 kl 



7P 4vr kl - v 2 
and for 77 collisions 



vbK x {yb) - -^-j [K {ub) - K (k o b)} ] (9) 



4f (6) = ^klbK x {bk ) (10) 

As for P^ ad , this is clearly expected to be 0(a em ) and from VMD one would expect 
1/250. It should be noticed that the eikonalised minijet cross-sections do not depend on 



A 77 and P 77 d separately, but depend only on the ratio of the two [^TJ j22|; which for our 



ansatze for A 77 means ratio of ko and P 77 d . From phenomenological considerations |l8|and 
fits to HERA data, fixing k Q = 0.735 GeV one finds a value P^ ad ss 1/200, which indicates 
at these energies a non-VMD component of ~ 20%. 

As mentioned, the QCD description requires the definition of ptmin- From HERA data, 
one notices that while lower values of Ptmin, i-e. 1.4 GeV, can be invoked to describe the 
beginning of the rise, a higher value, i.e. 2.0 or even 2.5 GeV, is better suited to describe 
the rise at higher energy. In Fig. la we show the fit to HERA data obtained with the above 
parameters, using a purely phenomenological fit to determine the non-perturbative part of 
the cross-section. 

Having thus established the range of variability of the quantities involved in the cal- 
culation of total photonic cross sections, we now proceed to calculate and compare with 
existing data the eikonalized minijet cross-section for 77 collisions. For photon photon 
collisions, we use the central value ptmin = 2.0 GeV. We also use P^ ad = 1/204, and A(b) 
from eq. (|10D with 3 different values of k Q which correspond to values within two standard 
deviations from the ZEUS [2y] collaboration value. Our predictions for 77 collisions are 



shown in Fig. (lb). A comparison with existing data shows that data points are better 
fitted by a higher larger value of ko, and we choose ko — 1 GeV. 

As stressed, the theoretical description is rather unsatisfactory and we now move to 
present a model for the overlap function which, in principle, should allow for a clearer 
predictability and to provide an expression for A (b) which could be applied to various 
cases of interest. We shall use Bloch-Nordsieck techniques to sum soft gluon transverse 
momentum distributions to all orders and compare our results with both the intrinsic 
transverse momentum approach as well as the form factor approach. In what follows, we 
shall first illustrate the Bloch-Nordsieck result and show that it gives a gaussian fall-off 
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Figure 1: Total inelastic photon-proton and photon-photon cross-section as described in 
the text. 



with an intrinsic transverse size consistent with MonteCarlo models |12] . We then calculate 
the relevant distributions and discuss their phenomenological application. 

In ref.([|23[]) it has been proposed that in hadron-hadron collisions, the b-distribution 
of partons in the colliding hadrons is the Fourier transform of the transverse momentum 
distribution resulting from soft gluon radiation emitted by quarks as the hadron breaks up 
because of the collision. This distribution is obtained by summing soft gluons to all orders, 



with a technique amply discussed in the literature |24], |25| . The resulting expression [pq , |27 
is 



Fw*(K±) = ^ / bdbJ (bK ± )e- h ^ M ^ (11) 



with 



7T 



9r* rM Hh, k 2 M + a/M 2 - k 2 . 



where cp = 4/3 and the hadronic scale M accounts for the maximum energy allowed in a 
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single (k 2 = 0) gluon emission. 

The definition given in eqs.(||]||), requires for its consistency a normalized b-distribution, 

i.e. 



d 2 bA(b) 



(13) 



so that the proposed Bloch-Nordsieck expression for the overlap function A(b), satisfying 
the above normalization, reads 



.4 



-h(b;M,A) 



BN 



2tt / bdbe- h ^ M ^) 



(14) 



An inspection of eq. (|T2"|) , immediately poses the problem of extending the known asymp- 
totic freedom expression for a s to the very small k± region. To avoid the small k± divergence 

and freeze a 



in eq.fll2l), it has been customary to introduce a lower cut-off in k 
i.e. to put 

oc s (k±) 



at k 



12tt 



1 



0. 



(15) 



with a = 2 in ref. 28 



33 - 2N f lnpi + a 2 A 2 )/A 2 )] 
For applications where the scale M is large (e.g., W-transverse mo 



mentum distribution calculations) eq.fll2l) is dominated by the (asymptotic) logarithmic 
behaviour and the small fc^-limit, albeit theoretically crucial, is not very relevant phe- 
nomeno logically. However, this is not case in the present context, where we are dealing 
with soft gluon emission in \ow-p t physics (responsible for large cross-sections). The typ- 
ical scale of such peripheral interactions, is that of the hadronic masses, i.e. we expect 
M ~ 0(1 2 GeV) and the small k± limit plays a basic role. This can be appreciated on 
a qualitative basis, by considering the limit bM << 1 of eq. (|T2"D. In this region, one can 
approximate 1 — Jo(kb) « b 2 k 2 /A, to obtain 



h(b; M, A) « b 2 A 



with 



A 



2L 

47T 



dk a. 



, k 
'A 2 ' 



In 



4M 2 



k 2 



(16) 



(17) 



One obtains a function h(b; M, A) with a gaussian fall-off as in models where A(b) is the 
Fourier transform of an intrinsic transverse momentum distribution of partons, i.e. 
exp(— k\/4:A 2 ). Note that the relevance of an integral similar to the one in eq.(|l7|) has 
been recently discussed in connection to hadronic event shapes [^9| . 

Our choice for the infrared behaviour of a s for a quantitative description of the distri- 
bution in eq. (|T2|) , does not follow eq. (|I5D , but is inspired by the Richardson potential for 
quarkonium bound states [ 30| . In a number of related applications [31, 32), it has been 



proposed to calculate the above integral using the following expression for a s 



a s (k A 



12tt 



P 



(33-2JV»ln[l+p(V)*] 
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which coincides with the usual one-loop expression for large (relative to A) values of k±, 
while going to a singular limit for small k±. For the special case p — 1 such an a s coincides 
with one used in the Richardson potential |30|], and which incorporates - in a compact 
expression - the high-momentum limit demanded by asymptotic freedom as well as linear 
quark confinement in the static limit. In [|31|] we have generalized Richardson's ansatz to 
values of p < 1. For 1/2 < p < 1, this corresponds to a confining potential rising less than 
linearly with the interquark distance r. The range p ^ 1 has an important advantage, i.e., it 



allows the integration in eq.(12) to converge for all values of k±. For the motivations given 
in [|3lj the value p = 5/6 was chosen in previous calculations of the transverse momentum 



distribution of Drell Yan pairs [31, 33 



Having set up our formalism, we shall now examine its implications. The distribution 
A{b) depends upon the hadronic scale M in the function h{b). This scale depends upon the 
energy of the specific subprocess and, through this, upon the hadron scattering energy. It 
plays a crucial role, just as it did for the Drell- Yan process, where the expression of eq. (|TT|) 
has been successfully |34|, [35], [33], ^] used to describe the transverse momentum distribution 
of the time-like virtual photon or W-boson. In these cases[[JJ], [531, the scale M was found 
to be energy dependent and to vary between y/CP/4 and \[Q 1 /2. In the calculation of the 
transverse momentum distribution of a lepton pair produced in quark- ant iquark annihi- 



lation [34 1, the scale M was obtained as the maximum transverse momentum allowed by 



kinematics to a single gluon emitted by the initial qq pair of cm. energy y/§ in the process 

qq^g + i(Q 2 ) (19) 

In the Drell- Yan case, one needed h(b) to calculate the transverse momentum distri- 
bution of the lepton pair, here we use it to evaluate the average number of partons in 
the overlap region of two colliding hadrons. In this case e~ h ^ is the F-transform of the 
transverse momentum distribution induced by initial state radiation in the process 

qq^ jet jet + X (20) 

where the jet pair in process (^) is the one produced through gluon-gluon or other 
parton-parton scattering with total jet-cross-section Oj e t and X can also include the quark- 
antiquark pair which continues undetected after emission of the gluon pair which stimulated 
the initial state bremsstrahlung. We work in a no-recoil approximation, where the trans- 
verse momentum of the jet pair is balanced by the emitted soft gluons. Then the maximum 
transverse momentum allowed to a single gluon is given by 

QmaAs) = ^(1 - § ~f) (21) 

with yj §j e t being the jet-jet invariant mass over which one needs to perform further inte- 
grations. An improved eq.([|) now reads 

n(b,s) =n soft (b,s) + Y1 J ~ J —fi( x i)fj( x 2) J dz J d Pt A BN(b,q max ) - (22) 
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where are the quark densities in the colliding hadrons, z = Sj et /(sxiX2), and -^jr^, is the 
differential cross-section for process (^0) for a given p t of the produced jets. 

Unlike the usual expressions for n(b, s), eq.(^) does not exhibit factorization between 
the longitudinal and transverse degrees of freedom since the distribution A BN depends 
upon the quark subenergies. Factorization can be obtained however, through an averaging 
process whereupon one can factorize the b-distribution in eq . (|22|) , by evaluating A BN with 
Qmax at its mean value, i.e. write 

n(b, s) = n soft (b, s) + A BN (b, < q max (s) >)a jet (23) 

with 



a 



jet 



^I^I^^hl^ (24) 



1,3 

and 



qmaX{S) " Ei,j I ^fi/aM J ^f j/b {x 2 ) J(dz) ^ 

with the lower limit of integration in the variable z given by z min = 4pj min /(sxix 2 ). To 
grasp the energy dependence of this scale, one can use a simple toy model, in which the 
valence quark densities are approximated by 1/ ^fx and thus obtain 

3 \fs 

< qmax(s) >~ -Ptmijn- (26) 

for 2p t min « \fs- For p t min = 1-4 GeV, as in typical eikonal mini-jet models for proton- 
proton scattering 0, one obtains values of < q ma x(s) > which range from 0.5 to 5 GeV 
for y/s between 10 GeV and 14 TeV respectively. A more precise evaluation of the above 
quantities depends upon the type of parton densities one uses, and will be discussed in a 
forthcoming paper. 

From the discussion about the large b-behaviour of the function h(b), we then expect 
ABN(b, s) to fall at large b more rapidly as the energy increases from y/s = 10 GeV into 
the TeV region. In Fig. 2, we compare this behaviour for the function A(b) with the one 
obtained through the Fourier transform of the squared e.m. form factor of the proton, 
eq.( ||). The function A(b) from the Bloch-Nordsieck model is calculated for A = 0.1 GeV 



and values of < q max > which include those obtainable from eq.(|26|) in the energy range 
y/s w 10 GeV -T- 14 TeV. 

We notice that, as the energy increases, A(b) from the form factor model remains 
substantially higher at large b than in the Bloch-Nordsieck case. As a result, for the same 
a j e t the Bloch-Nordsieck model will give smaller n(b, s) at large b than the form factor 
model and a softening effect of the total eikonal mini-jet cross-sections can be expected. 



In conclusion, we have studied the parameter dependence of the eikonalized mini-jet 
model for photonic and hadronic total inelastic cross-section and found that a large un- 
certainty arises through the description of parton distributions in impact parameter space. 
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Figure 2: Comparison between the A(b) distribution function from the Bloch-Nordsieck 
model (full) and the form factor model (dots). 

A model, derived from soft gluon summation techniques, is described and compared with 
expectations from the currently used form factor models for such disatribution. Such 
comparison indicates a distinctly different behaviour in this large b-region suggesting a 
softening of the rise of the total cross-section in mini-jet models relative to the ones with 
the hadron form-factors. 
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HCMP CT92-0026, C.S.I.R. (India), grant No. 03 (0745)/94-EMR-II and US-Department 
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